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(57) Abstract: The invention is directed to a method for diagnosing the state of a system. The system may be mechanical, chemical, 
electrical, medical, industrial, business operations, manufacturing related, and/or processing related, among other. The method may 
measure a signal from the system. Further, the method may compare the signal to an expected signal. The method may calculate a 
signal strength and/or a noise. The signal strength and noise may be functions of a frequency. Further, the signal strength and noise 
may be used to determine a channel capacity and/or a rate of information. A comparison ci the raie of information and the channel 
capacity may yield information associated with the state of the system. The information may be used in diagnosing the state of the 
system. Further, the expected signal may be derived from a model. The model may be tuned to the measured signal. The model may 
have parameters that are associated with features and/or faults of the system. These parameters may be used in diagnosing the stale 
of the system. Further, selectively repeated diagnosis over time may yield a prognosis of the system. 
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MODEL-BASED MACHINE DIAGNOSTICS AND PROGNOSTICS USING THEORY OF 
NOISE AND COMMUNICATIONS 

5 BACKGROUND OF THE INVENTION 

Related Applications 

This application claims priority of U.S. patent Application, Serial No.60/235,251, filed 
September 25, 2001 entitled: "MODEL-BASED MACHINE DIAGNOSTICS AND PROGNOSTICS 
USING THEORY OF NOISE AND COMMUNICATIONS", and is incorporated herein by reference in 
10 its entirety. 

1. Field of the Invention 

The present invention generally relates to a method of diagnosing systems. In particular, the 
present invention relates to a model-based method for diagnosing the operational health of a system, 
15 and for forecasting the future operational health of the system 

2. Description of Prior Art 

Diagnostics and prognostics are used in many fields. These fields may include mechanical, 
chemical, electrical, medical, manufactining, processing, and business operations, among others. Each 
of these fields has problems and difficulties relating to determining the source of a problem, identifying 
20 the severity of the problem, and predicting the behavior of a system in relation to the problem. 

For example, reliability and maintenance of complex equipment is critical to productivity and 
product quality. The purchase price of many typical equipments may account for half the equipments' 
cost Maintenance and support during the "life-cycle" may consume an amount roughly equal to the 
book value of the asset 

25 For example, at a typical chip plant billions of dollars are invested in equipment; many 

traditional manufacturing plants invest hundreds of millions of dollars in manufacturing machinery. 
Most maintenance is rigidly scheduled by "time in service", not condition. Machine productivity is lost 
during maintenance downtime. This and unscheduled downtime due to failures represent a very large 
part of a machine user's cost of operation. 

30 In these embodiments, machines are complex systems of components: gears, shafts, bearings, 

motors, lead screws, sensors, electronics, microprocessors, etc. integrated into a working whole. 
Machines in this context can also be biological, chemical, or hydraulic, among others. Defective or 
degraded components, alone or interacting, can render a machine dysfunctional. The rriachine may fail 
catastrophically and not complete its task, or it may lose tolerance, resulting in defective parts. 

35 Although methods and models exist for many individual component failures, errors from slightly 
degraded components can "stack", yielding overall system failure even when these models predict 
health of all individual components. 
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Many typical designers of machinery and engineers that maintain machinery, essentially know 
what the potential system faults are, and at what locations in the machine these faults will occur. 
However, unexpected breakdowns will still occur. Many typical designers and engineers do not know 
and/or can only poorly predict when fruits will occur. Further, periodically, healthy machinery must be 
5 taken out of service for maintenance. Thus adding an unnecessary cost Also, it may be very difficult 
and sometimes impossible to observe many of the conditions internal to the machine that lead to failure. 

One can easily imagine metaphoric extensions of these problems to other fields such as 
chemical, electrical, medical, rnanufactiiring, processing, and business operations, among others. As 
such, many typical systems suffer from deficiencies in providing accurate diagnostics and prognostics. 
10 Many other problems and disadvantages of the prior art will become apparent to one skilled in the an 
after comparing such prior art with the present invention as described herein. 
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SUMMARY OF THE INVENTION 

Aspects of the invention are found in a method to diagnose operational health of a system, and 
to forecast future health. For example, the method may permit intelligent scheduling of mamteuance 
downtime in a mechanical or chemical system. Further, the method may be used, for example, to avoid 
5 functional and catastrophic failures. 

Further aspect of the invention may be found in the method assembling models of the rrachme 
system, including system components and known system faults. Faults may be treated as "noise". In 
addition, parameters in me model may be "tuned" nrom signals from the real system, causing the model 
to mimic the real system in its present condition. Diagnosis may then be performed by observing the 
10 model. 

In another aspect of the invention, the method may treat the system as a communications 
c hann el estimate signal and noise levels, and diagnose health of the system with a tuned model by 
assessing how much information per unit time the system in its present condition can convey over its 
communications channel. The method may compare this maximum amount to the amount required by 

15 the system to execute a certain task. The method may assess if a system, in a given state, can perform a 
certain function, with a specified performance, within a given tolerance. 

Other aspects of the invention may be found in a method based on fundamental principles of 
physics and information theory. Further aspects may be found in the method assessing functional 
condition or state operable to perform a specified task, in addition to potential for catastrophic failure. 

20 Additionally, the method may operate on a tuned model, avoiding interpretation of complicated signals. 
Furthermore, the method may allow predictive scenarios for a system's possible future health and 
functional condition, given certain observed trends. 

In another aspect, for a different system or a new system design, only the model may be altered 
and not the basic diagnostic algorithm. The method may alsc permit the incorporation of knowledge of 

25 faults, and the intent of the designers of systems, into the diagnostics routines. 

Another aspect of the invention maybe found in assembling detailed dynamic systems models 
of the system in question. The models may posses a one to one correspondence between at least 2 
portion of the physical components or elements in the real physical system, and elements in the dynamic 
systems model. In some embodiments, all of the physical components are modeled in a one to one 

30 correspondence. 

In one embodiment the model may include all possible faults and potential failures in the 
system models. These effects may be tabulated as "noise" in the system. Noise in a signal is the 
difference between the actual signal and the expected signal. In the model, noise may be induced by 
changes in parameters of dynamic system elements, which then alters any signals passing through a 
35 system Alternatively, if a certain fault cannot be described by these means, then sources of noise may 
be inserted into the system model, at locations in the model that are consistent with the locations of the 
faults in the real system. The intensity of these noise sources may be adjusted to make the model behave 
like the real system 
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In an exemplary embodiment, the method may further include, placing sensors on a machine, 
to moniior the machine; exciting the machine; and observing the machine's response via the sensor 
outputs. The collected data may be used to rune the models parameters, so the model niimics the real 
system. After data has been collected on the actual rnachine, the system model may be excited with the 
5 same excitation as the actual rnachine. The outputs of the model may be compared to the corresponding 
outputs of the real machine. If the model's outputs differ from the real machine's outputs, the values of 
model parameters may be adjusted or changed, including the intensity of the noise sources, until the 
model's outputs approximate the actual system's outputs. 

In one ernbodiment, the channel capacity, C, of the system may be calculated. The channel 

1 0 capacity may be the maximum amount of information per unit time that can be measured from 
successfully conveyed through the machine. The channel capacity may depend on the design and 
construction of the system, and the present condition of the system, which results from aspects. These 
aspects may include manufacturing, aging and damage, among others. For example, faults maybe 
encoded as "noise" in the model. AnalyticaUy, the channel capacity may depend on the strength of the 

1 5 noise levels in the system, relative to the strength of the excitation system response signal. 

In the exemplary ernbodiment, for a desired task to be performed by the niachine, the rate of 
information associated with the task may be calculated. The rateof information may depend on the 
desired speed at which the machine does me task, the desired loads imposed on the machine, the 
complexity of the task, and the desired accuracy at which the machine should do the task. Further, the 

20 rate of iiuormation may be measured. 

Another aspect of the invention may be found in comparison of the rate of information to the 
channel capacity. This comparison may be used to evaluate the operabiliry of the system. If the rate of 
irftirmation is less than or equal to the channel capacity, the system may perform the desired task within 
the desired precision. If the rate of information is greater than the channel capacity, the system may 

25 functionally fail. 

Another aspect of the invention may be found in die formulation of extremely uemiled models 
of the system to describe a system's condition. In one exemplary embodiment of a system, the model 
includes bond graph based models of a motor, a gear box, and other mechanical transmission 
components. These extremely detailed models (a) exhibit a one to one correspondence between 

30 elements in the model and components in the real system; (b) incorporate many typical effects of the 
device into the model including defects; (c) include in the models, via finite element concepts instilled 
into bond graphs, the dynamically distributed nature of components in the real system, and (d) use noise 
sources to account for defects and degradation of components. Simulation of the motor and gear box 
models may generate the complex spectra measured during operation of these devices. These models 

35 may iruraic real system behavior and may be used to store information regarding the health condition of 
the machine. 
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In a further aspect, the models tabulate the effects of system faults (system maladies) as 
"noise" in the machine. Noise may be the difference between the actual signal received, and the 
expected signal that should be received As a machine degrades or ages, the difference between actual 
and expected signals may become larger. Thus noise levels may increase. These noise methods permit 
5 incorporation of faults into me models that heretofore could not be described analytically. The herein 
described methods have imported this body of knowledge to mechanical, hydraulic, other physical 
systems, and others, to name a few. 

In an additionally aspect, the method maybe used to predict the fixture conditions of systems, 
for scheduling maintenance and avoiding functional and catastrophic failures of the systems. The 
10 method may forecast if a complex system is capable of performing a given task, at a given speed and 
load, within a specified tolerance. 

The model system may be implemented on a computer system Hardware and software 
components may in combination allow the execution of computer programs associated with the method. 
The computer programs may be implemented in software, hardware, or a combination of software and 
15 hardware. 

Further modifications and alternative ernbodiments of various aspects of the invention will be 
apparent to those skilled in the art in view of this description. Accordingly, this description is to be 
construed as illustrative only and is for the purpose of teaching those skilled 

As such, a method for diagnosing and prognosticating the state of a system is described. Other 
20 aspects, advantages and novel features of the present invention will become apparent from the detailed 
description of the invention when considered in conjunction with the accompanying drawings. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

Figure 1 is a schematic block diagram depicting the Shannon-Weaver Model for use according 
to the invention. 

Figure 2 is a schematic block diagram depicting the information path according to the invention. 
5 Figure 3 is a schematic block diagram depicting a series of information paths according to the 

invention 

Figure 4 is a block schematic diagram depicting a computation system for implementing the 
method, according to the invention. 

Figure 5 is a schematic block diagram depicting a network system for implementation of the 
1 0 method, according to the invention. 

Figure 6 is a block flow diagram depicting an exemplary method according to the invention- 
Figure 7 is a cross sectional view of squirrel cage induction motor. 

Figure 8 depicts Ghosh and Bhadra's [5] bond graph of a squirrel case induction motor. 

Figure 9 depicts the stator resistances in Figure 8 redistributed to each of the stator coils. 
1 5 Figure 10'depicts a simplified representation of the signal and modulated GY element 

Figure 1 1 depicts a squirrel cage rotor with five bars. 

Figure 12 depicts a transformation of a and (3 phase currents into rotor bar currents. 
Fieure 13 depicts the bond graph structure including stator and rotor bar action. 
Figure 14 depicts the bond graph equivalence used in modeling. 

.«? — . 

20 Figure 1 5 depicts the bond graph representing stator and rotor bar action in the magnetic circuit. 

Figure 16 depicts angular velocity of rotor axis and stator currents in stator winding. 

Figure 17 depicts angular velocity of rotor axis and stator currents in stator windings, at startup. 

Figure 1 8 depicts angular velocity of rotor axis and 5-currcnts in each rotor bar, at startup. 

Figure 19 depicts angular velocity of rotor axis and 5-currents in each rotor bar, at startup. 

25 Figure 20 depicts angular velocity of rotor axis and 5-currcnts in each rotor bar, from startup to 

steady state. 

Figure 21 depicts stator currents and rotor velocity of a machine with a broken rotor bar. 
Figure 22 depicts stator current of 2nd phase and rotor velocity of a healthy rmchine at steady 

state. 
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Figure 23 depicts stator current of 2nd phase and rotor velocity of a machine with a broken rotor 
bar at steady state. 

Figure 24 depicts the angular velocity of rotor axis and 5 currents in each rotor bar when the 3*d 
bar is broken. 

5 Figure 25 depicts a torque-time plot of healthy rmchme and one rotor bar-broken machine. 

Figure 26 depicts rotor velocities of healthy and shorted machines. 
Figure 27 depicts rotor torques of healthy and shorted machines. 
Figure 28 depicts rotor bar currents of shorted machine. 

Figure 29 depicts Kim and Bryant's bond graph of an induction motor with state variables. 
1 0 Figure 30 depicts angular position and velocity of rotor axis- 

Figure 3 1 depicts flux in rotor <Z windings ; the P winding flux is similar. 
Figure 32 depicts flux in stator & windings ; the P winding flux is similar. 
Figure 33 depicts rotor velocity of a machine with a broken rotor bar. 
Figure 34 depicts stator current in the Frequency domain with broken bars. 

1 5 Figure 35 depicts torque-speed characteristics of the ideal and degraded machines. 

Figure 36 depicts power spectrum of the machine response and noise. 

Figure 37 depicts noise in the signal of the angular velocity of the degraded machine. 

Figure 38 depicts channel capacities with a broken bar. 

Figure 39 depicts rotor velocity of ideal and shorted machines. 
20 Figure 40 depicts power spectrum of angular velocity for the shorted machine. 

Figure 4 1 depicts spectral content of stator current of phase A; (a) Ideal machine, (b) Shorted 
rmchine. (c) Ideal machine of [15]. (d) Shorted machine of [15] 

Figure 42 depicts spectral content of stator current of phase A with two severely shorted coils. 

(Rsi = &2 =0.0079 a n s2 =10) 
25 Figure 43 depicts channel capacities with one shorted coiL 

Figure 44 depicts channel capacities with two shorted coils. 
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DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT 

Claude Shannon formulated a mathematical theory of communication. His groundbreaking 
approach introduced a simple abstraction, the conmumication channel consisting of a sender (a source 
of information), a transmission medium (with noise and distortion), and a receiver (whose goal is to 
5 reconstruct the sendefs messages), see figure 1 . 

The transmitter injects messages from an iriforrnarion source into the channeL The receiver accepts 
a signal from the channel thai contains the transmitted signal altered by the dynamics of the channel, 
and corrupted by noise added by the channel. 

An analogy is made between a machine component or a system and a cornmuni cations channel. 

10 During operation, information is sent as a signal over a communications channel from transmitter to 
receiver (See Figure 2). The signal over the channel is altered by limited dynamic bandwidth, 
nonlinearities and noise. The goal is for the receiver to extract and reproduce the message, dsspite 
distortions and noise. Design of communications systems is aided by powerful theorems of Shannon 
(1949), which establish minimum signal to noise ratios for error free transmission. 

1 5 A machine component (or system) accepts a "signal" from an upstream component, by its function 

alters that signal, and then passes the "signal" on to the next downstream component In the analogy of 
this article, a rnachine is a communications channel. When operating properly, the "signal" from an 
upstream component is "received" by a downstream component. Faults in the machine that disrupt 
functionality alter the "signal". Faults will be viewed as agents that alter system parameters or 

20 contaminate the signal with "noise". Unless the signal to noise ratio is kept sufficiently high, 
downstream components cannot "resolve" the "signal message" error free, and the machine 
malfunctions. 

In performing a nvection, a machine, . component or system accepts a stimulus "signal" from 
another upstream component, alters that signal via its mechanical- faction, and then passes the signal 
25 on to the next component. The signal contains information, which can be envisioned as a "message" to 
other components in the machine. The "message" relates to the function or intended operation of the 
rnachine or machine components. The mechanical ftmciion often includes kinematics of motion and 
dynamics of operation. 

Here we will strike an analogy between a machine component and a communications c hann el. The 
30 transmitter, an upstream component, activates our rriachine component "channel" with input signal x(t). 
Visage of the "signal" through the channel is associated with component functionality: component 
kinematics and dynamics alter the signal. The component response defines the outpui y(t). When the 
component channel operates properly, the "message" contained in the signal "received" by downstream 
components can be unambiguously "resolved" . Defects and degradation of the component afflict normal 
35 operation, "arstorting" the signal and contanainating it with "noise" n(t). Unlike electromagnetic 
communication channels, the signal may pass through multiple power domains: electrical, mechanical, 
solid, fluid, chemical, biological, etc. along its path through a component or a rnachine system. We can 
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view a machine as a chann el consisting of individual component channels connected together to form a 
larger channel. 

The theory depicted and described in this application may be adapted in unique ways not 
contemplated by others, including, but not limited to, Shannon's theory, A component is designed to 

5 have functionality, which can be defined in terms of the (designer's) intended reaction of the component 
to an excitation. Degradation alters the component response. For example, fretting corrosion of the 
surface of an electrical contact changes the electrical impedance through the contact Although this 
alters the response to a voltage stimulus, the resulting signal distortions caused by changes of the 
electric contact impedance are often posed in terms of an effective noise riding on the transmitted 

10 signal. Thus degradation of the contact via fretting is often modeled as an effective noise source and/or 
an impedance change. 

In communications theory, Shannon's theorems traditionally estimate the maximum rate of 
information C that can be transmitted through a communications channel, given its bandwidth w and 
ratio of signal to noise powers S/N. Designers of traditional conimunicanons channels considered C to 

15 be fixed, and their designefforts focused on designing transmission or encoding schemes that would 
increase the rate of information R up to its upper limit, the channel capacity C. If applied in a 
nontraditional manner to machinery, Shannon's theorems can yield a threshold signal to noise ratio 
(S/N)t, In the communications channel analogy, dynamics inherent in the component functionality can 
be included in bandwidth w and in the signal to noise ratio S/N of the channel capacity C. These 

20 dynamics may change as the component degrades, causing C to change. For a typical machine and 
components, signal transmission rateR should be constant, since rnachme or component operation is 
often repetitive (or penodic) and at or near steady state: the machine controller and/or upstream 
components continue to inject their signals into a machine (or component), regardless of its condition. 
When applying communications theory to a machine component or system, we will first trace the 

25 path of signal power flowing through a healthy (functional) component or machine, to define the 
communications channel through the component or machine. Along the signal path we will list the 
various forms of energy oi H o*ver into which the signal is transformed. Functionality will be defined in 
terms of the input to output response, for components or the system: if for a given set of input 
excitations, the output response matches within some tolerance the desired output, the component or the 

30 system is functional; otherwise it is dysfunctional. If needed, we will consider each separate 
energy/power domain and its transduction as a communications channel, and then connect these 
channels together in a manner consistent with the machine's functionality and design. Bond graphs 
(Kamopp, Margolis, and Rosenberg, 1990), which map power flows through dynamic systems, can be 
useful, since bond graphs readily handle systems with diverse energy domains in an energetically 

35 consistent manner. 

After analyzing the healthy system, we will then incorporate component faults and degradation 
modes into the system model To affect functionality, the degradation effects must alter or block the 
flow of signal power through the component. Questions we must answer include: How does each 
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degradation mode alter the signal flow, and affect system or component parameters? Does the particular 

degradation cause components to become nonlinear? Does the particular degradation generate another 
signal, Le., noise? We will incorporate degradation into the system model as changes to existing system 
model parameters or as additional elements (e.g., sources of noise). Location of each degradation mode 
5 in the system model will be consistent with the locus of the degradation in the physical machine 
component oi system. 

Aspects of this method include: 1) Individual components, or an entire machine system consisting 
of multiple components, can be analyzed. 2) System malfunctions can be predicted, including individual 
faults and those due to a collection of seemingly healthy components. Errors from slightly degraded but 
10 individually healthy components can stack through a machine system, rendering it unable to meet 
tolerance. 3) The current status of the system, and time to system malfunction can be estimated by 
simulations based on these models. 

Other aspects of the diagnostic procedure include: Deterrnine and trace the path of the signal flow 
through the healthy system, from signal in to.signal out. For the sick system, model the faults with noise 
15 sources or parameter changes. Multiple system outputs may exist At each output, rally the signal power 
and the total noise power to obtain a signal to noise ratio S/N. Estimate the bandwidth w for the signal 
path through the degraded component communications channel, using the enhanced system model. 
Apply Shannon's theorems to diagnose the absolute health of the machine component or machine 
system The health of each individual component in a machine syrtem can be assessed, and likewise the 
20 health of the entire machine system. 

The analysis of each machine component communications channel may contain the following: 
. Haaithy Machine Model which has no raults and functions perfectly. This is an idealization that 
reflects the machine designer's original intended concept. The output y 0 (t) of a signal x(t) 
propagating through this ideal machine will define the intended message or signal y 0 (t) mat the 
25 machine component or machine system communications channel is supposed to transmit and 

receive. The signal powers S - P{y}and Si « P{y 0 } defined by S - P{y}= 7 [ D>(0f*«pi (1), 
both in the channel capacity C and the rate of information transmission R through the channel, will 
be based on this ideal machine. The resulting model is simple, the concept of perfect health is well 
defined, and the signal y 0 (t) that the receiver is supposed to receive is well defined. 

30 

. MachiDfi EaultS : These include common degradation raults for a given component Common 
examples include pitting of gear teeth, fatigue cracking of shafts, and deterioration of insulation on 
electric motor stator or rotor coils. 

35 . Machine Eailli Models incorporate the Machine Faults as sources of noise n(t) and/or changes 
in system parameters consistent with impenections, faults and degradation modes of a particular 
machine element. Noise will be defined as any signal component that should not be in the perfectly 
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transmitted and received message signal y 0 (t). This may includes harmonics generated by nonlinear 
elements. 

. n^flrarfpri Marhin** Model : This is the overall system model that results from adding the 
5 Machine Fault Models to the Healthy Machine Model It includes sources of noise n(t) and changes 

in system parameters. When all noise sources are zero, the healthy machine results. Transmission 
of the signal x(t) through the degraded machine (noisy cornmunications channel) induces received 
signal y*(t), generated by signal x(t) (sent through as y 0 (t))and noise n(t). 

10 The analogy may also be extended to a set of machines, a process, a manufacturing or assembly 

method, or others. The analogy may hold for a series of "information channels" as seen in figure 3. 

The model-based diagnostics is based on fundamental first principles of physics and 
information theory. The methods uses sensor signals to tune the parameters of a model of the system, 
such that the model men mimics the operation of the real system. Diagnostics are performed on the 

15 model. The diagnostic system can be designed as part of the design of a new machine. Also, models 

allow what if predictive scenarios for a machine's possible future health and functional condition/ given 
certain observed trends in the machine's health. For a different machine or a new design, only the 
model of the operation of the system must be altered, not the basic diagnostic algorithm Models also 
avoids interpretation of complex sensor signals, trying to figure out what a particular peak or dip, or a 

20 band of frequencies in a signal means, in terms of machine health. Instead, time wise changes to 
machine parameters can be followed, and projection of these trends can be used to forecast future 
health. Models also permit incorporation of knowledge of faults, and the intent of the designers of 
machinery, into the diagnostics routines. 

To quantitatively analyze transmission through the channel, Shannon introduced a measure of 

25 the amount of infonnation in a message. The measure is related to the probability of occurrence of the 
events for which the messages are about A message that informs the receiver that a rarely occurring 
event is about iu happen contains the most information. A message iriforming about an ^ready "known" 
event conveys little information. Information entropy, a measure of the average amount of information 
(or uncertainty) in a message, can be denned [ 1] as 



30 H = \ 



~ X Pi *° S 2 Pi tf x * s discrete 

0) 

- p(x) log 2 p (x)dx if x is continuous 



Here Pi is the probability of occurrence of the message's event xi if the random variable is discrete and 
p(x)is the probability density function for the random variable if the random variable x is 
continuous. Here Pi is the probability of occurrence of the message's event x, xi if the random variable 
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is discrete, and />(*)is the probability density function for the random variable *, if the random 

variable .t is continuous. 

Shannon's entropy rate (R) measured a source's information production rate, and the channel 
capacity (Q measured the information carrying capacity of the channel. As per one of Shannon's 
theorems [1], if ^ ^ C, then there exists a coding technique which enables transmission over the 
channel with an arbitrarily small frequency of errors. This restriction holds even with bounds the noise 
in the channel. A converse to this theorem states that it is not possible to transmit messages without 
errors if R > C . Thus the channel capacity is defined as the rnaximum rate of reliable information 
transmission through the channel. 

In another theorem, Shannon derived the channel capacity for a time continuous channel with 
additive white Gaussian noise. His expression 



10 



C = a>log 2 



(2) 

involves the average transmitter power, 

1 0 

1 5 of me signals * c (0 , the power of the noise, 

N = P{*t)}=±)[n(t)]dt (4) 
and bandwidth CO of the channel in hertz. If the bandwidth is non-fiat, then the capacity of the channel 



is given by 

V. < 5) 



c=[.„ 8 ,|, + .W 



M' w>, 

20 Similarly, the entropy or inforniation rate for messages 

R -^log^S,/^) (6) 

derived by Shannon involves S { , the average power of the desired signal to be transmitted, A,- , the 
maximum aUowed RMS error between recovered and original messages, and <0 { , the signal bandwidth. 
Shannon's communication theory could be applied to tne fault diagnosis of machine systems. 

25 A machine component (or system) accepts a signal &om an upstream component by its function alters 
that signal, and then passes the signal on to the next downstream component. In Bryant's analogy, a 
machine conveys infonnarion in a signal and is thus a communications channel. When operating 
properly, the signal passes through the system and is successfully received within desired tolerances at 
the machine's output. Faults that disrupt operanon alter the flow of signal. Paults will be viewed as 

30 agents that contaminate the machine's signal with "noise". Unless the signal to noise ratio (S / 2f) is 
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kept sufficiently high, downstream components cannot resolve the signal message error free, and the 
machine malfunctions. 

Noise is defined as an "unwanted signal tending to ohscure or interfere with a desired signal", 
as "any signal which interferes with the transmission of a signal through a network or tends to mask the 
5 desired signal at the output terminals of the network", and as "an unwanted signal tending to interfere 
with a required signal". Thus noise is the difference between the actual signal received, and the signal 
desired to be received. To apply this definition to mechanical systems, we must define the desired 
signal. We shall call this desired signal the "ideal" signal x o (t) . Note that x o (0 , an idealization, must 
be produced by a system without noise. This is possible only with models, not with real systems. 
10 The "ideal" and "degraded" models may be defined as follows; 

• The ideal machine model has no faults and functions perfectly. Its output defines the signal * a (0 that 
the machine channel is supposed to receive. From this, we can estimate signal power S i . 

• The degraded rnachine model is the overall system model that results from adding faults to the model. 
We will incorporate faults as noise n(t) . Thus the signal *(0 = x o (t) + n(t) contains noise 

15 defined as any signal component that should not be in the perfectly received message signal. Noise is 
any deviation from the ideal signal, including unwanted harmcnics generated by nonlinear elements. 
This will estimate the noise power N . 

We can incorporate degradation or imperfections into the system model. Degradation can be instilled 
in a bond graph model by varying bond graph parameters, adding noise (effort or flow) sources, or 
20 changing the power pathways. 

The models may take various forms. These forms may be any form appropriate for use in the system 
"of application. For example, these forms may be heuristic, neural networks^detenninistic, probabilistic, 
and others. 

The method and model system may be implemented on a computer system, S see Figure 4). 

25 The term "computer system" as used herein generally describes the hardware and software components 
that in combination allow the execution of computer programs. The computer programs may be 
implemented in software, hardware, or a combination of software and hardware. A computer system's 
hardware generally includes a processor, memory media, and input/output (I/O) devices. As used 
herein, the term "processor" generally describes the logic circuitry that responds to and processes the 

30 basic instructions that operate a computer system The term "memory medium" includes an installation 
medium, e.g., a CD-ROM, floppy disks; a volatile computer system memory such as DRAM, SRAM, 
EDO RAM, Rambus RAM, etc.; or a non-volatile memory such as optical storage or a magnetic 
medium, e.g., a hard drive. The term "memory" is used synonymously with "memory medium" herein. 
The memory medium may comprise other types of memory or combinations thereof. In addition, the 

35 memory medium may be located in a first computer in which the programs are executed, or may be 
located in a second computer that connects to the first computer over a network. In the latter instance, 
the second computer provides the program instructions to the first compuicr for execution. In addition, 

013270.00015:124601.01 



WO 02/27418 



PCT/US01/29935 



~ R 
the computer system may take various forms, including a personal computer system, mainframe 

computer system, workstation, network appliance, Internet appliance, personal digital assistant (PDA), 
television system or other device. In general, the term "computer system" can be broadly defined to 
encompass any device having a processor that executes instructions from a memory medium. 
5 The memory medium preferably stores a software program or programs for the reception. 

storage, analysis, and transmittal of information produced by an Analyte Detection Device (ADD). The 
software program(s) may be implemented in any of various ways, including procedure-based 
techniques, component-based techniques, and/or object-oriented techniques, among others. For 
example, the software program may be implemented using ActiveX controls, O+ objects, 7avaBeans ; 

10 Microsoft Foundation Classes (MFC), or other technologies or methodologies, as desired. A CPU, such 
as the host CPU, for executing code and data from the memory medium includes a means for creating 
and executing the software program or programs according to the methods, flowcharts, and/or block 
diagrams described below. 

A computer system's software generally includes at least one operating system such Windows 

1 5 NT available from Microsoft Corporation, a specialized software program that manages and provides 
services to other software programs on the computer system. Software may also include one or more 
programs to perform various tasks on the computer system and various forms of data to be used by the 
operating system or other programs on the computer system The data may include but is not limited to 
databases, text files, and graphics fries. A computer system's software generally is stored in non-volatile 

20 memory or on an installation medium A program may be copied into a volatile memory when running 
on the computer system. Daia may be read into volatile memory as the data is required by a program. 

Further, the method may be implemented across a set of networked devices (See Figure 5). 
The method may be performed remotely from the system Further, the results of the method may be 
transmitted, stored, processed, and accessed across a neiwork, among others. 

25 For example, parameters for a model of a patient's health may be stored on a smart card. 

These may be accessed anHj-n.mbined with the method to deterrnine a charge in state of the patient's 
health. In another exemplary embodiment, a machine may be located in a remote location. A service 
provider may periodically access data from the machine from a remote location and diagnose the 
machine. These diagnoses may be used in predicting the failure of the machine. Further, these 

30 diagnoses may be used in placing an order for a replacement. 

Figure 6 depicts a flowchart for diagnosing according to the invention. The method may be 
implemented in software and/or hardware. Further the method may include some or all of the steps in 
various combinations. 

In a first step, the user is directed to assemble detailed dynamic systems models of the machine 
35 system in question. The models may possess a one to one correspondence between physical components 
or elements in the real physical system, and elements in the dynamic systems model. One may include 



01 3270.00015: 124631. 01 



WO 02/27418 



PCT/USO 1/29935 



15 

all possible faults and potential failures in the system models . This invention may tabulate the effects of 
faults as "noise" in the system. Noise in a signal is the difference between the actual signal and the 
expected signal. In the model, noise may be induced by changes in parameters of dynamic system 
elements, which then alters any signals passing through a system Or, if a certain fault cannot be 
5 described by these means, then sources of noise (often white noise) will be inserted into the system 
model, at locations in the model that are consistent with the locations of the faults in the real machine. 
The intensity of these noise sources can then be adjusted to make the model behave like the real 
machine. 

One may then judiciously monitor the machine or system Excite the machine or system, and 
10 observe the machine's or system's response, for example, via the sensor outputs. 

One may then tune the model' s parameters, so the model mimics the Teal system Excite the 
system model with the same excitation as the previous list item Compare the outputs of the model to 
the corresponding outputs of the real rmchine or system If the model's outputs differ from the real 
machine's or system's outputs, adjust or change values of model parameters, including the intensity of 
15 the noise sources, until the model's outputs closely match the actual system's outputs. 

One may then manipulate the model, which now rntmics the real machine or system in its 
present condition: 

From the model, one may calculate the channel capacity, C. of the machine. C is the luaximum 
amount of information that can be observed successfully conveyed through the machine. The channel 
20 capacity depends on the design and construction of the system, and the present condition of the system, 
which results from manufacture, aging and damage. Faults are encoded as "noise" in the model. 
Analytically, C depends on the strength of the noise levels in the system relative to the strength of the 
excitation system response signal. 

For a desired job to be performed by the machine, one may calculate the rate of information R 
25 associated with the job F depends of the desired speed at which the machine does the job, the desired 
loads : the complexity of the job, and the desired accuracy at which the rnachine should do the job. R is 
measured in bits of information per second. 

- Compare R to C. If R < C, the machine will perform the desired job within the desired 
precision. If not, the system has functionally failed. 

30 The comparison of R to C may yield a diagnosis. Alternately parameters of the tuned model 

may yield a diagnosis. Further, this diagnosis may be associated with the determined noise. In 
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addition, the noise and/or diagnosis may be indicative of combined faults. Further, combined variances 
in parts, while within tolerance limits, may comprise a fault, defect, or others. 

The method may be repeated over time to build a prognosis of the machine or system For 
example, a prognosis may predict the failure of a part. 
5 Further, the method may be applied to many systems such as those depicted above. In 

addition, parameters from the tuned model may indicate the type or state of a defect, fault illness, or 
condition, among others. 

In typical applications, the method may involve formulation of extremely detailed models of 
machine devices to describe a machine's condition. These are critical to success. For example,, included 

1 0 are bond graph based models of a motor, a gear box, and other mechanical transmission components. 
These extremely detailed models: (a) exhibit a one to one correspondence between elements in the 
model and components in the real system; (b) incorporate all known effects of the device into the 
model, including defects; (c) include in the models via finite element concepts instilled into bond 
graphs the dynamically distributed nature of components in the real system, and (d) use noise sources to 

15 account for defects and degradation of components. Simulation of the motor and gear box models can 
generate the complex spectra measured during operation of these devices. 

The models tabulate the effects of system faults (machine maladies) in a very novel way: as 
"noise" in the machine. Noise is the difference between the actual signal received, and the expected 
siraal that should be received. As a machine degrades or ages, the difference between actual and 

20 expected signals becomes larger, and thus noise levels increase. These noise methods permit 

incorporation of faults into the models that heretofore could not be described analytically. The concept 
of noise has been used he a vily in electronics and communications engineering, to design around noise 
''faults 1 ' always present in these electronic and electromagnetic systems. Electronic noise, including 
resistor noise, shot noise, burst noise., and flicker noise among others has been generally tabulated or 

25 modeled with noise sources placed in a model of the electronic circuit This work imported this body of 
knowledge to mechanical, hydraulic, and other physical systems, but in addition, systems extended the 
modeling schemes of noise to include noise induced by changes in parameters of the system. 

The method also applies techniques of information theory to rjiachinery - as opposed to present 
applications that are limited to electronic communications systems - to quantitatively assess the current 

30 health state of a machine. The method treats a machine, such as a CNC engine lathe, as a noisy 

communications channel, to assess reliability and functional condition. A message transmitted and 
received over a cornmunications channel picks up noise due to imperfections present in the physical 
channel. For example, music transmitted over an AM channel is overwhelmed by buzzing when the 
receiver is near electrical power transmission lines: the transmitted musical message is obscured at the 

35 receiver by electrical noise. In an analogous manner, a machine transmits a message over a machine 
channel. For example, a lathe, viewed as a communications system, has transmitter = CNC controller, 
channel = (drive motor + gear box + lead screw, + tool carriage on ways + cutting tool / workpiece). 
and receiver = workpiece "Noise" includes effects of fatigue, spurious vibration (from other rmchines), 
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. and other errors due to wear of machine and cutting components. A transmitted "message" is properly 

"received" if the finished part is within tolerance, or in a general rnachine, if the machine performs its 
function within specified tolerances. Excessive noise in the rnachine system may cause a part to be out 
of tolerance, or causes the general machine to operate outside the specified tolerance limits. With this 
5 view, Shannon's cornmumcations theorems may be applied to machinery-. Shannon's theorems may 
accurately estimate the limits on the amount of information per unit time C that can be sent through a 
noisy communications channel. C depends on the channel's state, including dynamics and signal to 
noise strengths (ratios). For a lathe, making a part of certain geometric complexity at a given speed, to 
within a desired (fidelity) tolerance is characterized by an information rate R. IfR. < C, Shannon's 

10 theorems predict success; if R > C, the part will be out of tolerance. As a machine deteriorates, C 
decreases, and eventually R > C. Now the machine cannot make the part with the same speed and 
tolerance. The channel concept appears to be a very sensitive discriminator of a machine system, even 
for the stacked effects of a collection of moderately degraded components. 

The method may be used for predicting the future conditions of machinery, for scheduling 

1 5 rnaintenance and avoiding functional and catastrophic failures of said machinery. The method can 

forecast if a complex system is capable of doing a given task, wimin a specified tolerance. A multitude 
of parameters associated with the machine's model may be tuned, such mat the model emulates the real 
system. 

These modeling and system assessment techniques could be useful to designers of rmchinery, 

20 to assess the efficacy, reliability and durability of a design under various user conditions. 

In addition to mechanical systems, these methods could apply to almost any kind of dynamic 
system, including chemical, electrical, medical, manufacturing and processing, and business operations, 
among others. For example, in the medical world, a detailed model describing the dynamics of the 
caitiio- vascular system could be developed. This model would possess multiple parameters that 

25 describe behavior and condition of the heart and blood vessels, and their interactions with other body 
systems such as lungs and kidneys. The model could be tuned from medical signals and data derived 
from tests and procedures, such as Electro-Cardiogram, blood pressure, and data from lau tests and 
radiology. After tuning the models, a channel capacity C could be estimated to assess the condition of 
that system, and compared to a rate of information R This comparison would assess the health state of 

30 the patient. The rate of information would describe the ability of the cardio-vascular system to perform 
at various levels characterized by tasK speed, load, complexity, and tolerance. Since the rate contains 
these factors, degrees of health and sickness could be assessed quantitatively or assessed, in a formal 
manner. This could automate medical diagnostics. Medical prognostics would extrapolate trends of 
parameters in the model, or trends contained in the data, and apply the channel capacity and rate of 

35 information concepts of communications theory, to forecast future health scenarios. 

These methodologies could be extended to evaluate business practices, procedures, and 
enterprise structures. A business operation has dynamics imposed by its processes, people, and 
structure. The application would treat an enterprise as an imperfect corroiunications channel, and 
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construct models of information flow through that system. Transmitters— the orders— will send 
information over imperfect "enterprise communication channels". Imperfections— problems in the 
enterprise, or interference between coiiflicting missions— adds "noise" to channels. Receivers— the 
customers— must receive the message— the product— within tolerances— customer expectations— 
5 despite noise. The application would define "channels" through enterprise units, construct models that 
mimic these channels, and then apply communications theory to diagnose and prognose these channels. 

The models in these embodiments and claims can take various forms: from structured methods 
such as bond graphs, differential equations, and finite elements, among others, to heuristic methods 
such as neural networks, fuzzy logic, expert systems, and other computer methods. 

10 The method for applying communication theory to machines and systems need not be limited 

to signals derived from models. The method could be extended to signals measured from real systems. 
Here the ideal signal Xo(t) could be approximated from measurements taken from a real machine, or 
from several rnachines, in excellent condition. The difference between x c (t) and the signal x(t) measured 
from a degraded machine could replace those derived from models, mentioned earlier. Similarly. ±e 

15 difference could be used to confirm that a noachine operates within tolerances. Further, an ideal signal 
could be a signal from a machine with a known defect. The difference between the signals would then 
confirm a specific defect, among others. 



Fx^mr^ry application tn a squirrel canfi induction motor 

20 Equation numbers in this example refer to equations listed in this subsection. Similarly, an 

appendix is attached mat is referenced in this subsection. 

One exemplary application of the invention is a method for diagnosing an induction motor. 
For example, a motor has two major sub systems: a rotating rotor and a static staror. Induction machines 
can have a wound rotor, or a squirrel cage rotor. Widely used squirrel cage induction rmchines exhibit 

25 great utility for variable speed systems and are simple, rugged, and inexpensive. The squirrel cage rotor 
is a structure of steel core larninarions mounted on a shaft, with solid bars of conducting material in the 
rotor slots, end rings, and usually a fan. In large machines, the rotor bars may be of copper alloy, dri vcn 
into the slots and brazed to the end rings. Rotors of up to 50 cm diameter usually have die-cast 
aluminum bars. The core larninations for such rotors are stacked in a mold, which is then filled with 

30 molten aliiminum. In this single economical process, the rotor bars, end rings and cooling fan blades are 
cast at the same time. 

Figure 7 is a schematic of a squirrel cage induction motor. A substantial literature modeling 
induction motors employs Park's (1929) two-reaction theory, which accounts for magneto-mechanical 
energy transduction via multi-port inductances. From Park's model, Ghosh and Bhadra (1993) 
35 formulated the bond graph in Figure 8. We altered Ghosh and Bhadra's bond graph to partition and 
make explicit the electrical, magnetic, and mechanical energy domains; to form a one to one 
correspondence between physical components in the machine, and elements in the bond graph; and to 
append additional elements to the bond graph to make it more consistent with real induction motors. 
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When energized by an AC supply voltage, the stator coils form a radial magnetic field vector 
that rotates within the interior of the stator, about its central axis. Within this interior the stator field 
cuts through the squirrel cage rotor, including conductor bars that extend axially. This time varying 
field induces a voltage over the rotor bars. Resulting bar currents flow in the sequence: bar end ring 
-» opposite side bar -» opposite end ring h> original bar. Induced by this time varying current loop is a 
secondary magnetic field, which attempts to align with the stator field. However, because the rotating 
stator field induced the secondary field of the rotor, the stator field leads the rotor field, and 
consequently, the rotor chases the stator field, always following. This is motor action (Lawrie, 1987). 
The induction motor speed depends on the speed of the rotating stator field. 

The real system we will consider is a two pole, 'V connected three phase squirrel cage 
induction motor. In (Ghosh and Bhadra, 1993; Sahm, 1979; and Hancock, 1974), a muiti phase 
induction motor was modeled with an equivalent two-axis representation. Each phase winding generates 
its own magnetic field, which can be represented as a vector aligned along the axis of the winding. The 
sum of these phase vectors produces a phasor vector. If the phase vectors vary properly with time, the 
phasor rotates. 

A transformation from three phases (a,b,c) to two phases (a,p) was represented in (Hancock, 
19-74) in matrix form. If the V and 'a' phase windings are co-axial, the induced Magneto Motive 
Forces (MMF) of the V and 'a' phases of the three and two phase systems are co-directionaL By 
appropriate changes to the two phase currents, the magnitude of the phasors of me three and two phase 
systems can be made equal. Ghosh and Bhadra (1993) represented this in their bond graph via 
transformer elements in the stator section. The two phase currents were represented in terms of three 
phases as 

"V2/V3 -i/v6 -l/ys] ' 

0 l/V2 -1/^2 j ' 
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"cosO cos2.t/3 cos4,t/3 
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>- 
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(1) 



Under assumpnons of a spatially sinusoidal distribution of MMFs, and ignoring magnetic losses and 
saturation, Ghosh and Bhadra (1993) expressed a symmetric mducuoh motor in an orthogonal 
stationary reference frame with a and p phases fixed on the stator as 
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Equation (2) relates stator voltages 'to stator and rotor currents. In addition, needed is the electro- 
magnetic motor torque for a P-pole machine, expressed as 

Pr. 



This motor torque is balanced against other torques via 



' 2 



(3) 
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(4) 



Terms on the right side of equation (4) represent rotor inertial torque, shaft/bearing damping 
torque, and load torque, respectively. In equations (2) to (4), V aS and V p , are a and p axis stator 
voltages; i^ and ip s are a and p axis stator currents; iar and i Br arc a and 3 axis rotor currents; R s and R r 
5 are stator and rotor resistances; L,, and L, are stator self inductance, mutual inductance and rotor self 
inductance; T e and Tl are electro-magnetic torque and mechanical load torque; J is the moment of 
inertia of the rotor } c is the viscous resistance coefficient; cor and aim are electrical and mechanical 
angular velocities of the rotor; and P is number of pole pairs. 

Ghosh and Bhadra (1993) represented equations (1) to (4) in their bond graph, reproduced in 
1 0 Figure 8. They used modulated gyrators MGYrr^Lmip,, MGY.T 2 =Lripr, MGY:r 3 =L ra i a3: MGYT4=L r i cl r, 
to represent the electro-magnetic torque of equation (3); employed transformers TF:mt, TFmfc, TF:nk, 

TF:m4 ; TF:m 5 with moduli ■ m, = J^, m 2 = 7tz 3 = -^6 , m 4 = -Jl, m 5 = —Jl to implement the 

mathematical transform of equation (1); and excited the system with effort sources MS e :v a , MS e rvb : and 
MS c :v c having sinusoidal voltages with equal amplitudes but 0, jc/3, and 2r.B phase lags, respectively. 
1 5 Although this correctly programs the governing equations for a three phase induction motor, it lacks a 
correspondence between bond graph elements and real system components. Moreover, elements and 
their constitutive laws involve only electrical and mechanical energy do m ains. Faults or design 
parameters relevant tc the magnetic domain are only implicit in the mutual inductances, posed as 2 port 
inertances I:a and I: (3 with constitutive laws 
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In ?. C 5> A-ps, Xao and A^ r are flux linkage of the respective windings. In Figure S. five integral 
(independent) causalities exist on inertance energy storage elements, with system state variables Xas, 
Xps, *. a r» Xp r> and h, where h is the rotor angular momentum 

To represent real system elements or components explicidy, certain bond graph elements 

25 should be moved, altered or added. In Figure 8, a and (3 phase stator resistance elements, Rsa and R s0 
should be split into three stator coil resistances R 3a , R*, and R 3C) without alterating the governing 
equations. The revised bond graph shown in Figure 9 moved and R^ back through the transformers 
in front of the phases. To maintain an equivalence between Figure 8 and Figure 9, we must relate R sa , 
R t b and Rsc to R^ and R i{J . Since most motors possess symmetry between phases, let R M - R& = R S c = R, 

30 and R Sa - R 50 = R s . For the bond graphs of Figure 3 and Figure 9 to be equivalent the voltages (efforts) 
to the 2-port inertances on the stator sides must be equal in both Figure 8 and Figure 9. The causality in 
both Figure 8 and Figure 9 asserts that these voltages to the 2-port inertances arise from the neighboring 
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]- junctions. Summing voltages from other bonds to these 1 -junctions, and equating these respective 
voltages between Figure 8 and Figure 9 gives 



m A m s m A [ 

m 2 m 3 



yjn 2 m A )\ m s 



in, m 




r 



W2 2 



\V b -R 



K™2 n h 4 
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5 By solving for ip*/^, we obtain equations in terms of resistances and transformer moduli 



(8) 



_ mfm;m;m A m 5 R s - (rrqni^m^ + m?n£m A m s + m{n£m 4 m 5 )R 
(mfm 2 mlm 5 + m^m^m.m 4 )R 

By replacing the transformer moduli, m~m 3 of the three phase to two phase transformation with real 
numbers, m, = m 2 = m 2 = S, m 4 = 4l, m. = -V 2 , which is given in equation (1), we find that 
R 5 =R, i.e., R Sa = R sp = R„ - R s b = R 5 c 

Simplified Representation of the Signal and Modulated GY Elements 

In terms of the 2-port I field of equation (5), equation (3) can be rewritten as 



(9) 



From this relation, Figure 8 can be rearranged into the form of Figure 10 : where the modulated gyrators 
MGY: in = \p 2nd MGY: n 2 = Aor are modulated by the flux linkages ^ and ^ r of the 2-port 
inertances. 

The number of squirrel cage rotor bars depends on the rotor's size, and usually, tens of bars 
are in one rotor. In this study we consider the squirrel cage rotor with five bars (numbered 1 to 5) 
depicted in Figure i 1 . Shown also is the rotor magnetic field (dashed line), with north poles (N) on top 
of the rotor, and south poles (S) beneath, and bar currents. Currents directed out of plane are denoted by 
a V, and currents flowing into the plane are denoted by a V. Each end of each rotor bar is attached to 
a solid end ring. Induced currents flow through each bar and end rings. With five bars, there exist five 
different currents (flows) in this rotor. At the instant of the rotor position shown in figure 1 l-(a), the 
sums of the currents induced by the rotating magnetic field of the stator in bar 1, 2 and 5 must be equal 
to the sum of the currents in bar 3 and 4. Likewise, the current surnrnation of bar 1 and 5 at the position 
of Figure 1 1 -(b) must equal the sum of currents h bars 2, 3 and 4. In Figure 1 1 , the thickness of each x 
and • shows the relative current magnitude in each bar. 
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To incorporate individual rotor bars into the bond graph, the a and p phase currents and 
voltages of the rotor should be split into separate bar currents and voltages. The a, b, c and a, p axes 
are stationary with respect to the stator, but because the rotor rotates relative to these axes, bar currents 
must depend on the rotation angle 6 of the rotor. Using results in Hancock (1974), rotor bar currents 
can be related to the a, p phase currents as 

In equation (10), i rk represents the current in the k* rotor bar (k = 1, 2, ... n), X* and fy t are rotor 
currents from Figure 8, and magnitude modulus m depends on the total number of bars, n. For n - 5 
bars, we will have currents i r! to h<. Accordingly: rotor bars can be incorporated into the bond graph of 
Figure 10 via 0-modulated transformers. 

Figure 12 shows the transformation of a and p phase currents into individual rotor bai 
currents, where the transf onner moduli are 



mr k = mzos<0 + 



mr V4 



= msinjfl-!- f 



k = l f 2 f ...,n 



with n = 5 . 



(11) 



(12) 



(13) 



In Figure 12 : the battery of 0-junctions on the right side completes the summation of a and p phase 
currents demanded by the right side of equation (10). The voltages that sum over the two 1 junctions 
located between the I fields and the MTFs give rise to 

K + (™M + (™2 + • • ■ + (mr s = 0 
-V + ('"O^; + (mr-)^ + • ■ ■ + O^ioHo = 0 
Here the flux linkage \ u X 2j ■ ■ A.10 associated with rotor bars are located to the right of the MTF's. To 
obtain the torque contributed by each bar, equation (10) for k - 1, 2, ., 5 is rewritten in matrix form as 
cos 6 sin 9 

cos(0 + -y) sin(0 + ^) 

cos<0 + — ) sin(0-r— ) 
= m 5 5 

COS< # -r — ) Sini 0 t — ) 

The two column vectors of the 5x2 transformation matrix A form an orthogonal set for any value of 
rotor rotation angle 6; the rank of A is 2. For the mxD (mDn) matrix A having rank n, there exists 
(Strang, 198S) an nxm left-inverse B such that BA=I D , where I Q is the identity matrix of order n. In our 
model 



=Ai, 



'two phase 



(14) 



A'A = ;» i - 



a 5 ri 0' 



2 L° V 



(15) 



013270.00015:!24601.01 



WO 02/27418 



PCT/US01/29935 



23 

2 2 

and the left-inverse of A is A T if m = -j , i.e., the transformer modulus m has a value which 

normalizes A T A. For a rotor of n bars, m = J~ . The proof is shown m Appendix for this subsection. 
From equations (14) and (15), the inverse transformation is 



41 



cos9 cos(0 + -^) cos(0 + ^) cos(&A costfA 

J J «2 J 

sm0 sin(^-r^) sin(04-y) sin(0+y) sin(^ + -y) 



(16) 



10 



15 



If substituted into the rotor output torque equation (9), the electromagnetic torque becomes 



(17) 



(18) 



The revised bond graph in Figure 13 includes stator and rotor bar interactions based on equaiion (17). 
Here the moduli of the k* modulated gyrator is 

2(fc-l>^ 
n >J 

where n = 5 for Figure 13. Finally, the electric resistances of the rotor were grouped with each iotor bar 
in £ manner similar to that of the stator resistances. 

The bond graph in Figure 13 models the interaction between stator coils and rotor bars with 2- 
port I elements — inductances— in the electrical energy domain. An inductance only describes storage of 
magnetic energy. Neglected are power losses and leakage effects in the magnetic domain, which may be 
caused by component deterioration. To describe these interactions, we replace all I inductance elements 
with equivalent combinations of gyrators and C elements, without violating causality. Figure shows 
equivalent bond graph representations between an"l"and a GY and C combination; and a TF and GY 
combination. 



20 In Figure 14, n is the gyrator modulus (the effective number of coil rums); m is the transformer 

modulus; I is the flux linkage; <j) is the magnetic flux [Wb]; M is the magneto motive force [A]; p is 
the permeance of the magnetic circuit element [H]; ci and a are efforts; and fi and ii arc flows. In 
Figure 14 -(a), through the gyrator relations /i = n(p and ni = M . Using the constitutive law of the C 
element, M ■ §lfj, the two port I elements pertaining to the a and p phases were converted into 2-pon 

25 C elements that now represent interactions between magnetic flux and magnetomotive force of the 
stator and rotor. Figure 15 shows the new bond graph with five rotor bars and the GY - C - GY 
combination that replaced the 2-port I. The gyrators were then moved through the bond graph to new 
locations more consistent vvidi motor components. The GY to the left of the 2-port C was moved into 
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the electrical section, where it now represents the action and number of rums of the stator coils. The GY 
leap-frogged the transformers that were based on equation (1), changing moduli of these transformers 
according to Figure 14 -(b). The GY to the right of the 2-port C skipped over a 1-junction, converting 
that 1-junction into the 0-junction shown in Figure 15. Similarly, a 0- and 1-junction to the left of the 2- 
port C in Figure 13 were converted to a 1- and 0-junction in Figure 15. In the bond graph of Figure 15, 
electrical energy inputs, transformation of energy from electrical domain to magnetic domain, 
mathematical phase transformations, power interactions between stator and rotor bars in terms of 
magnetic flux and magneto motive force, and mechanical rotor output are all represented and labeled 
In Figure 15, the two sets of gyrator moduli ris and n r stand for the effective coil turns which relate 
electrical and magnetic variables of stator and rotor, respectively. 

State equations were derived from the bond graph of Figure 15 with n s i = n$2 = n*3 = Q;, Rm - 
Rib == Rsc " R*. In terms of magnetic variables, the state equations i 



> are 



n s m \ n * m i n s n h 
n s m A n s m 5 

h 

<Par = -Rrn s 2 K-P„<Pas + P.<P.) + Vp, —J 

k 



(19) 



m 6 J 



m 6 ( Pi P r -P n r — ™ 
15 wh?re the magnetic state variables are stator and rotor phase fluxes <|> a <. <j>p5.- <jw» (h pr . and rotor angular 
momentum h. The constitutive law of the 2-port C element is 



(20) 







'p. 


5] 




JP°r. 




A 




\m w . 


'<?/*' 




P* 


"m £ 


**. 




A 


Pr. 





In the state equations, 
1 



44 -4 2 n s \ 2 P,P r -P m 2 



(21) 



20 where the permeances P. = -~T > P m ~ — >JL ->P r ~ — 2* are expressed in terms of coil turns and 

n; n s n r n r 

inductances of stator and rotor. Here n s is the number of effective stator coil turns, n r the number of 
effective rotor coil turns, § the magnetic flux [Wb], M the magnetomotive force [A], P is the Permeance 
[H], and h the angular momentum [N-m-s=kg-m 2 -sec]. 
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Simulations of a squinel cage induction motor used the bond graph simulation tool, 20-SM 
(Control Lab Products, 1998). For integration of state equations, a Ruuge-Kutta 4 th order method was 
adopted. Values of the system parameters for the simulations are presented in 

Table 1, some were identical to those used by Ghosh and Bhadra. 

5 Table 1 System parameters of a two pole, three phase squirrel cage induction motor 



10 



R sa , R 5b; R SC [£XI 


Stator coil resistance 


0.07S8 


R:l>Rr2,...,Rrl0[Q] 


Rotor bar resistance 


0.0408 


v a , Vb, v c [V] 


Input voltage amplitude 


230 


[Hz] 


Input voltage frequency 


60 


L 5 [H] 


Stator inductance 


0.0153 


MH] 


Rotor inductance 


0.0159 


Lm[H] 


Mutual inductance 


0.0147 


n 


Number of rotoi bars 


5 




Number of effective stator coil turns 


100 


n r 


Number of effective rotor coil turns (bar) 


1 


c [N-s/'m] 


Mechanical resistance 


0.15 


J [kg.m 2 ] 


Mechanical inertia 


0.4 



Shown in Figure 16 and Figure 17 are piols of rotcr angular velocity and stator currents versus time. 
The rotor velocity rises slowly to a steady state value of about 377 rad/sec; the stator currents oscillate 
at the input frequency with initial large amplitude. After about 1.5 seconds, the motor reaches steady 
state: the currents in stator windings decrease to a steady value and no oscillation of rotor velocity 
exists. Figure 16 plots the rotor axis angular velocity vs time when 230V, 60Hz three phase AC 
voltages are input to the stator coils. Theoretically, when 60Hz alternating inputs are given to a two 
pole AC motor, the output velocity should be 360CRPM (377 rad/sec) and the simulation yields a 
steady state value very close to this (the difference is due to the mechanical resistance load). Figure 17 
15 expands the Figure 16 time scale to show the three stator currents with 120° phase difference, duiing 
motor start-up. 

Figure 18-20 shows the currents in the five rotor bars and the rotor velocity. Recall there exists 2n/5 
phase difference between currents in neighboring bars. This is clearly shown in Figure 19, which 
represents the motor starting moment. While the 60 Hz frequency of the stator currents generate a 
20 constant rotational velocity of the rotating magnetic field, the frequency of currents in the rotor bars 
decrease continuously as the rotor velocity increases. This is related to 'slip' in induction motors, the 
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normalized difference between the electrical angular velocity of the aii gap MMF established by the 
stator currents, and the electrical angular velocity of the rotor (Krause and Wasynczuk, 1989). Slip is 
defined as 



(22) 



5 where co s is the synchronous speed, or the speed of the stator currents, and o> is the speed of the rotor. 
The magnitude and frequency of the currents and voltages of the rotor depend on the relative velocity 
between the rotating magnetic field and the rotor. In these simulations, this relative velocity maximizes 
at t = 0, where the slip is unity. As the rotor velocity increases, the relative velocity and the slip 
decrease, suggesting that the decrease of amplitude and frequency of rotor bar currents in Figure 18 are 

10 probably due to the decrease of slip. If o s = slip s - 0 and no current is induced in the rotor bars 
(hence no torque). However, the. steady state currents of the rotor bars in Figure 18 are not zero, (even 
though there is no external load) because of the Actional load of the bearing modeled as a resistance 
R:c in Figure 1 5. If an external load is applied to the motor axis, the slip should increase and therefore 
the current and voltage in the rotor bars should also increase. Figure 20 shows the currents in the rotor 

15 bars during steady state. 

All simulation results shown above are for a healthy motor. When rotor bars break, currents, 
velocity, and torque will deviate. Because we have a one-to-one correspondence between bond graph 
elements and machine components, it is possible to represent broken rotor bars by increasing the rotor 
bar resistance R r . In modern squirrel cage induction motors, bars and end-rings contact the rotor core. 

20 Due to this available current shunt, currents in a broken bar are not zero (Manolas and Tegopolous, 
1997), i.e., the resistance is not infinity. Figure 21 shows the stator currents and rotor velocity for a 
rotor with the third rotor bar broken. During the transient rise time, the rotor velocity increases, and 
exhibits oscillations. Even at steady state, there exists periodic deviations of rotor velocity. -With these 
deviations, the amplitude of the currents in the stator coils also change. These changes are more clearly 

25 presented in Figure 23; for comparison, a corresponding healthy machine simulation is shown in Figure 
22. Figure 24 plots the currents in each rotor bar, with bar 3 assumed broken. From Figure 24, the 
induced currents are largest in the two rotor bars nearest the broken bar. Figure 25 compares the torque 
characteristics of the healthy rnachine and broken bar machine. The rotor torque oscillates in the broken 
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bar machine, even at steady state. During startup, the oscillation of torque is larger in the broken bar 

machine than the healthy machine. 

Simulations of an induction motor with a short circuited stator coil are shown m Figure 26-28. 

In these simulations, the resistance of the shorted coil decreases, and the coil current, the magnetic 

5 fields, and the induced currents in the rotoi bars also change. Figure 26 shows a difference in rise time 

of rotor velocity between the healthy machine and the stator coil short-circuited machine. Figure 27 

shows the rotor torque for both healthy and shorted machines. The overall trend of the torques are 

similar, but there exists small amplitude and relatively high frequency oscillations in the short-circuit 

case. These oscillations are also seen in the rotor bar currents, Figure 28, compared with the rotor bar 

10 currents of me healthy machine, shown in Figure 16. 

A bond graph model of a squirrel cage induction motor was constructed, based on a prior bond 
graph by Ghosh and Bhadra (1993), that exhibited a one-to-one correspondence between the bond 
graph elements and real system components. Included were stator coil windings for three phases, 
mathematical transformations to incorporate two reaction theory, magnetic state variables to represent 

1 5 magnetic interactions between stator and rotor, individual rotor bars and contributions to the total rotor 
torque and velocity, and mechanical inertias and resistances. The simulations in this article had five 
rotor bars. Using this model, simulations of a healthy machine were compared to simulatioos of 
rnachines with a broken rotor bar breakage and a shorted stator coil. The degraded rrachine simulations 
predicted oscillations in currents and angular velocities, seen in real motors. 

20 Most induction motor designs employ three phase excitation of the stator. For a rotor with 

more bars, the bond graph of Figure 15 can be easily altered. More rotor bars can be included in Figure 
15 by adding additional pairs of power pathways to the right of the 2-port C's, such that n power 
pathways fan out from both a and p rotor phases. For the new value of n, these power pathways must 
update equations (1 1) and (12) governing moduli mi\ for the modulated transformers MT£:mr k and 

25 equation (IS) governing modulus t\ of the modulated gyrators MGY:r*. To update the 
electromechanical torque, in equation (17) we must replace the 5 in the upper index of the sum and the 
square root argument in the denominator with the new value of n. 

A Q^nnH FvpmplaTy Application 
30 This subsection refers to equations 1-6 ir. the detailed description. In addition, the remaining 

equation numbers refer to equations within this subsection. Further, an appendix is attached which is 
referenced in this subsection. 

In a further embodiment of an induction motor, the bond graph model of a squirrel cage 
induction motor from above is adjusted. This model includes stator windings for 3 phases, two-reaction 
35 theory, magnetic interactions between stator and rotor, individual rotor bar contributions to rotor torque 
and velocity, mechanical inertias, and resistances and losses. Although this model does not include 
certain critical phenomena of the induction motor - s.g., magnetic field with rotor eccentricity or rotor 
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dynamics - this model is simple and can illustrate bow to apply Shannon's communication theory to 

machine systems. 

In the system shown in Figure 29, MSe:V, , MSe:V b and MSe:V c indicate the 2-phase 
alternating voltage applied to the motor. The resistor element RjI^ models resistive losses in the stator 
5 windings of the motor. The gyrator GY:n s models the transition from the electric to the magnetic 
domain of the power flow in the system The modulus of the gyrator n 5 equals the number of turns of the 
stator coil. The battery of transformers TF:mk convert the 3-phase into a rotating phasor vector. The 
two-port capacitance elements C represent the interaction between stator and rotor fields. 

In the rotor, electric voltage is induced in the metal bars by time varying flux cutting the bar 

1 0 circuits. This represented as the battery of gyrators, which have moduli n, related to the number of turns 
of the rotor. The modulated transformers MTF:mr k relate angular position of the rotor relative to the 
flux field. The resistor elements R:R, represent resistive losses in the rotor circuits. Modulated gyrators 
MGY:r k convert bar currents on the rotor bars into torque; this is the magneto-mechanical interaction. 
The moduli for these gyrators depend on fluxes in the rotor. The final transformer TFrmm is related to 

15 the number of magnetic poles in the system Power lost by bearing friction is accounted for by the 
resistance R:c. The remainder of the power drives the output shaft. 
Moduli in the bond graph (Figure 29) are given as follows [6]: 
1) Moduli of three phases are (mk) : 

/T ill 

20 2) Constitutive laws for two port C fields arc : 

N (?) 

where, i ■ = and 

*< = 7^-rr • *« = t^" ~h ' *' = rr~h : 

L s is stator self inductance, L M is mutual inductance and A. is rotor self inductance. The gyrator 
25 moduli is the number of stator coil turns, and gyrator moduli n r is the number of rotor coil rums. 
3) The modulated transformers MTF:mric are : 

mi\ =,i-cos^g+ [ k = l,2,...,n 



X" 













'In 



in [ 



9 + — — > with n = 5 



(9) 



where n is the total number of bars. 
30 4) Tae moduli of the modulated gyrators MGY:r k are : 

r k = 1l r [<Pfir mr k " 9cr mr t+* J * = U, (10) 
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5) The modulus for transformer TFrm,,, is : 



m m =— : P p is number of poles 



(11) 



State equations were derived from me bond graph (figure 29) with 7i 5l = n s2 = »,j - "* »n - »ri<r- n r . 
The results of state equations are 



1 h 



5 ^=-7 



(12) 



Par =-{ mr * R rl +™r 2 2R rl + W R r3 ^ ^ R r4 + R r5 )^~ 

-(ni^mr^ -r mr z mrJL n + mr z mr z R rl + mr 4 mr 9 R r4 + mr 5 mr w R r5 y^ 

x 1 * 

* A = -(" r 6^n + mr, 2 R r2 + m7- 8 : i? r3 + mr 2 R ri + mr l0 2 R rS )Lf 

71 r 

-(m ri mr 6 R n +mr 2 mr.R^ +mr, -mr s R r3 +mr A mr 9 R r4 +mr s -mr l0 R r J-^ 



(13) 



M. 



10 



- (mr 6 ■ r, + m?-, • r, + mr t ■ r 3 + mr 9 ■ r A + mr l0 • r 5 ) — - 



1 h 



n. -m,,, J 



(14) 



l ff, A, 



2 \ 



-jlf + — M, 



4 y; 



l- 



2 A 



77 '/L 



=' 



1 



m A -n s 



2 \ 



c 2 , D 2 



—M„+—M t 



f 1 



J ;j s 2 + J2 s; " V'": 



(15) 



K - 



72, 2 +J?. 



J3 V m 3 



772. 



(16) 



/i = (j\ ■ mi\ + r 2 • ;72/% 4 r 3 ■ mr 2 + r 4 • mr 4 + r 5 - 7?/r 5 ) ■ — M „ 

-rfc -mr 6 +r ? mr 7 + r 3 -mr 8 + r 4 -mr 9 +r 5 -mr lQ ) — M A -—c 



(17) 
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where the magnetic state variables are rotor angular position 0 o and momentum h and stator and rotor 

phase fluxes 9m > 9* » 9ar and 9 fir • 

Simulation of a squirrel cage induction motor employed MATLAB® ! s Runge-Kutta 4 th order 
method with a time step St = 1CT 5 seconds. Values of the system parameters presented in Table 2 

5 were given by Kim and Bryant [6, 7]. 

Using this model simulations were performed for an ideal machine, which has no faults and 
functions perfectly according to designer's specifications, and a degraded niachine. The ideal machine 
will serve as a reference of desired d\Tiamic behavior. The degraded motor will exhibit common 
degradation modes, including rotor bar breakage and stator coil shorts. We will excite the ideal and 
10 degraded machine models with identical test signals, record these signals., and then estimate the noise as 
the difference between degraded and ideal machine responses to the same test signal. 



Table 2 System parameters of a two pole, three phase squirrel cage induction motor. 



Parameter 


Description 


Ideal 


Degraded 


Rrj> Rib> R*c 


Stator coil resistance [0] 


0.07SS 


0.00079 - 0.0709 


Rrl, ... ,-RriG 


Rotor bar resistance [Cl] 


0.0408 


0.0412 - 4.0800 




Input voltage amplitude [V] 


230 




f 


Input voltage frequency [Hz] 


60 




*, 


Stator Reluctance [l/H] 


5.85xl0 6 




*, 


Rotor Reluctance [l/H] 


563 




*. 


Mutual Reluctance [l/H] 


5.41x10* 




n 


Number of rotor bars 


5 




Th 


Number of effective stator coil turns 


100 


1-90 


n r 


Number of effective rotor coil turns (bar) 






c 


Mechanical resistance [Hs/m] 


0.15 




J 


Mechanical inertia [kg-m 2 ] 


0.1 




m. (-P/2) 


Number of poles □ 2 


2 





15 

Figure 30~32 show sample simulation results for a nominal or ideal motor, i.e., a motor 
without faults. These simulations arose from the model of equations (7) to (17), with the "Ideal" 
parameter values of Table 2. Plotted are selected motor state variables versus time, begmning with 
motor startup, i.e. the motor voltages were switched "on". The rotor velocity rises slowly to a steady 
20 state value of about 377 rai'sec, as the momentum and all other state variables reach steady state. 
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Theoretically, when 60Hz alternating inputs are given to a two pole AC motor the output velocity 
should be close to 3600RPM (= 377 rad/sec). Rotor and stator fluxes of & phases are shown in figure 
3 1 and 32; the P fluxes are similar. Flux amplitudes increase to steady state, consistent with the 
angular velocity. 

5 Various faults can be developed in motors. For example, stator coil shorts cause overheating, 

increasing core losses [8]; rotor bar breaks or cracks in the die-cast rotors cause very large electrical 
resistance [6, 7, 9]; and bent or cracked shafts make- the rotation wobble [10]. 

In this article, we will focus on a broken rotor bar, and shorted stator coils. When rotor bars 
break, steady state velocity and torque of the rotor will deviate from the ideal response. With the bond 
10 graph shown in Figure 29, a broken bar can be incorporated into the model by increasing selected rotor 
bar resistances R,. The range of deviation is given in Table 2, last column Figure 33 (upper curve) 
shews the response of the motor with a broken bar after being switched "on". In this case, the rotor bar 
resistance was increased 10 times from its norninal value of 0.0408 □,ohms, to R r = 0.408 □.ohms. 
Piotied is the angular velocity versus time, from start up. Here the angular velocity increases, during a 
15 transient time characterized by deviations of rotor velocity. 

Figure 33 (bottom curve) shows the simulated startup (step) response for a motor with a rotor 
bar having resistance increased 100 times, to R r » 4.08 □.ohms. This curve shows increased and 
persistent oscillations, compared to Figure 30 for the ideal machine. 

It is weli established that when rotor faults occur, rotor harmonic fluxes are produced which 
20 induce currents in the stator at frequencies of f[k/(P p /2)-(l-s)± s\ m Here / is the supply 
frequency, P p is the number of poles , k =1,2,3 ,4. . and s is slip defined as [1 1, 12] 

(18) 

CO 

s 

In equation (IS) 0) 3 is the synchronous speed derived from the frequency of the stator currents, and 

CQ r is the angular speed of the rotor. Slip can have a value from 0 to 1 . 
25 Figure 34 (a), constructed by applying a Fourier transform to the steady state portion of the 

simulation results of Figure 33 (bottom curve), shows some frequencies of the stator current of phase A 

in the vicinity of the excitation frequency (6OH2). Figure 34 (b), scanned from reference [13], shows 

spectral densities of typical currents versus frequency measured from a motor with three broken bars. 

Comparison of Figure 34 (a) and (b) show similar shape and location of spectral peaks. 
30 Figure 35 compares the torque characteristic of the ideal rnachine and broken bar rrachine R r - 

0.408 □.ohms. The rotor torque oscillates whenever the rotor velocity oscillates. Due to the rotor 

asymmetry' the level of pukating torque is increased [6, 12]. 

The average power in a signal x{t) , of duration T can be estimated as [14] 

S = P$c(t)}=± : )\x(0]dt (19) 
1 0 
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or as 



x 2 P°i 



If X„ is a sequence sailed from x(t) at equally spaced discrete instants. The power spectral density, 
the magnitude squared of the Fourier transform of x(t) , is given by 

wh=re X(f)=)x(t)e- Ja 'dt . 

—a. 

For discrete x„, we employed a fast Fourier transform to obtain X k . 

The total power can be calculated in the frequency domain, or in the time domain by Parseval's 
theorem. [14] 

,0 )\ X (tfdt=)\X(jfdf (22) 

-O -00 

The discrete form of Parseval's theorem is defined as [14] 

iw 2 4£w 2 - • (23) 

Equation (5) can be rewritten as 
15 C = ]{log 2 (S + AO - log : N}dco . (24) 

0 

Combiranon of the ordinal signal power spectral density (S) and the noise power spectral density (A 1 ) 
represents the signal power spectral density 

S. = S + iV P5) 

from the Jegraded machine. 

20 Shannon [1] assumed a Gauss^n white noise statistically independent of the signal. To remove 

this restriction, we need to calculate the noise power directly from the time domain signals. In the time 
domain, the noise is defmed as the difference between actual and ideal signals 
n(t) = Mt)-x o (t). (26) 
Here Jt(0 is the output of the degraded niachine, and x 0 {t) is the output of the ideal niacbine. As 

25 demonstrated in the Appendix for this subsection, removal of the independence restrictions between 
signal and noise admits negative channel capacities. 

Power spectral densities SI and N r can be defined as the magnitude squared of the Fourier 
transforms for signal x(t) ano noise "(0 respectively. To calculate the channel capacity with these 
values, we must replace (S + N) m equation (24) with Si , and N with AT , to have 
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Figure 36 (upper line) shows the power spectrum of the rotor velocity of the ideal machine as 
shown in Figure 30 (upper line), and denned in section 4.1. This figure was constructed by applying a 
fasi Fourier transform to the angular velocity data of Figure 30 (upper line). In this case, we assume 
5 zero noise, and thus the system functions perfectly, according to the design specifications. Using 
equation (27), we get an infinite channel capacity for an ideal machine system, since by definition, the 
noise and noise power are zero. 

If there are faults such as broken rotor bars as mentioned in section 4.2earlier, the power 
spectrum will change as noise contarninates the signal. Figure 37 shows the startup response x(t) of 

10 the machine with a broken bar (upper line), and the noise in the time domain from the degraded 
machine, defined by equation (26). This noise (lower curve and magnified in Figure 37) is the 
difference between the degraded machine's response curve in Figure 37, and the startup response of the 
ideal machine in Figure 30. The presence of several frequencies is evident. Figure 36 (dots in the upper 
line) shows the power spectra (signals and noise) of a degraded inachine, with a cracked rotor bar. Note 

15 that the power spectrum of the degraded machine signal x(t) nearly overlaps the power spectrum of 
the ideal machine signal x o (t) ; in the figure, the two almost coincide. In our model, we increased the 
resistance of broken (cracked) bar by 10%. Using equation (27), we obtained a channel capacity of 
1.1 xlO 6 (bits per second). Here the integraticn bandwidth CO in equation (27) was equated to the 
entire sampling bandwidth (5xl0 4 /7z) based on the Nyquisfs sampling rate, where St = 10 " 

20 seconds was the time step employed in the numerical solution routine. In this procedure, we viewed the 
numerical solution's data points as a "sampled" signal, with sampling interval equal to the numerical 
method's time step. The Nyquist's sampling rats gives the smallest bandwidth associated with the 
sampling interval S t = 1 0 5 seconds. 

In equation (6) for entropy rate R, S, represents the average power of the output signal from 

25 the healthy rnachine and iV,. represents the largest acceptable deviation, i.e., a tolerance on the noise. 
The signal bandwidth (#,-) was equated to CO; see the previous paragraph for justification. The 
functional requirements of the rnachine deterrnine the noise or error tolerance N i demanded by the 
machine to work satisfactorily. For example, if we have an application wherein the maximum allowed 
error or tolerance must be within 10% of the signal of the ideal niachme, and if we employ the same 

30 bandwidth as for the channel capacity, then from the equations (6) and (19), the rnformation rate [R) is 

f |K A ] 



R = w i log 2 (5,./^.) = 5xl0 4 log : 



\T J J 



= 3.3xl0 5 (bits per second) (28) 
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With a channel capacity of 1.1x10° (bits per second), the result for R above satisfies the condition of 
R^C . This suggests a still "healthy" machine. 

If the resistance of the broken bar increases 21 times to Rr « 0.8568D, then the channel 
capacity drops to 2.4 xlO 5 (bits per second) below the required (R) of 3.3x1 0 5 (bits per second). 
5 Since this result doesn't satisfy R^C, according to Shannon's theorem, the machine is 
malfunctioning. As the magnitude of the fault (bar resistance) increases, the channel capacity 

diminis hes. 

Figures 38(a)-3S(d), show selected power spectral densities of the stator current of phase A at 
steady state, for selected bar resistances. These figures are similar to Figure 34 (a) and are often used as 

10 diagnostic indicators. Side bands are absent for smaller values of Rr, but start to appear after the rotor 
bar resistance equals approximately 0.7670Dohms (1780%). Figure 38(e) also plots the channel 
capacities versus the percent change of bar resistance from the ideal value gfven by table 2. The dashed 
line indicates the 10% noise power tolerance (R =3.3 x 1 0 3 ) estimated in the previous paragraph. Here 
percent change is defined as %R bar =(R b ar -JO /J? L , where R bar is the current value, and 

15 Rl,. the ideal. The channel capacity of point (d) in figure 38(e) has a negative value; see the Appendix 
for reasons. As shown in figures 38(c) and 38(d), significant side bands with large intensity begin to 
appear in the power spectra, wherever the channel capacity curve sinks below the 10% information rate 
line (dashed curve). From a practical standpoint, for industrial grade rugged machines such as motors, 
we would begin to notice errors when these exceed 10% or more in the motor's output velocity. Thus, 

20 10% was chosen as the critical velocity condition. 

The curve of figure 38(e) can be separated into regions with three distinct slopes: region 0, 
which connects the infmiTe channel capacity of the ideal system to that of "real" systems; region I, with 
stable C and "healthy" operation (region I would be associated with the normal life cycle operating 
region of the system); and region II, where C declines to the (dashed) failure line. Note that the marked 

25 change in the slope of C or the rapidly dirnirnshing values of C, going from region I to 1% could presage 
iaihire. Figures 38(bV(d) suggest that once side bands appear, the slope of C becomes noticeably more 
negative. 

Signal based diagnostic methods, that trigger upon detection of side bands, at earliest would 
notice the broken bar fault at point (b) in figure 38(e); figure 38(b) suggests that detection of the tiny 
30 side band would be difficult. In contrast, the channel capacity curve's knee - where the slope abruptly 
changes - occurs at 1500%, before 1780% of figure 38(b). Here the abrupt change in slope might be 
easier to detect 

Simulation of an induction motor with short circuited turns on its stator coil is shown in figure 
39. Here stator resistance F*i of phase A was decreased 50%, and the effective number of turns 
35 represented by gyrator modulus n s was similarly decreased from 100 to 50. In this simulation, only one 
of the stator coils has shorted turns. In the model and physically, as turns are short circuited, the 
resistance in that coil decreases, and the effective number of turns also decreases. Figure 40 shows the 
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various signal and noise power spectra of this shorted machine. Note again that the power spectra of the 
ideal and degraded machines nearly coincide. The shorted coil seems to affect the angular velocity 
relatively less than a broken bar. The two startup responses in the figure 39 are nearly the same, 
especially at steady state. Essential differences at steady state are in the signal's phase, generally not 
5 contained in power spectra. By taking the difference between degraded and ideal response, the 
information on phase difTerences is conveyed in the noise power spectrum, in addition to the magnitude 
information. 

When stator coil turns short out, we observe only a rise in some of the frequency components 
which already exist in the stator current spectra of an ideal irachine [15]. Figures 41 (a) and (b) shows 

1 0 spectral content of the steady state stator current of phase A, from simulations of the bond graph model. 
For comparison, spectra from Gojko and Penman's model [15] are also shown as figures 41 (c) and (d). 
Figure 42 shows spectral content of stator currents for two shorted coils, phases A and B. In figures 43 
and 44 are plotted the channel capacities versus percent change in the coil resistance, for shorting of 
phases A, and A and B, respectively. Again inforrcation rate for the 10% noise level on angular velocity 

1 5 is shown as the dashed line. 

Similar to the broken bar case of figure 38, figure 43 exhibits a ''healthy" region I. with stable 
channel capacity, and a region II with sharply dimi n ishing channel capacity. Again the sharply changed 
slope of region II could prognose failure. 

In this article, it was demonstrated how Shannon's theory of communication could be applied 

20 to machineiy, to utilize Shannon's powerful theorems. Concepts of rate of information rate and channel 
capacity were reviewed, and applied to an induction motor. At the heart of the method is machine 
4i noise 1 \ estimated as the difference between actual and ideal responses. By subtracting the ideal 
- response x o (t) from the acmal response *(0, the noise signal contains only rnformation about the 
faults. From the noise and signal were calculated power specUa, used in equation (27) for channel 

25 capacity. Rate R from equation (28) serves as a critical values dependent on the system's tolerance to 
errors, here called "noise". The channel capacity was calculated for a motor with shorted stator coils 
and broken bars. The channel concept agreed with other existing fault monitoring methods, but results 
suggest that it could detect faults much earlier. It can be concluded from this study that the channel 
capacity concept could serve as an effective Qjscrirninator of motor and machine system health. 

30 It is to be understood that the forms of the invention shown and described herein are to be 

taken as the presently preferred embodiments. Elements and materials may be substituted for those 
illustrated and described herein, parts and processes may be reversed, and certain features of the 
invention may be utilized independently, all as would be apparent to one, skilled in the art after having 
the benefit of this description of the invention. Changes may be made in the elements described herein 

35 without departing from the spirit and scope of the invention as described in the following claims. 

As such, a method for diagnosing the state of a system is described, in view of the 
above detailed description of the present invention and associated drawings, other 
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modifications and variations will now become apparent to those skilled in the art. It should 
also be apparent that such other modifications and variations may be effected without 
departing from the spirit end scope of the present invention as set forth in the claims which 
follow. 
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Nomenclature 

C channel capacity 

c viscous resistance coefficient 

/ supply frequency 

5 H information entropy 

h angular momentum [N-m-sec] 

J moment of inertia 

Ls, Lm, Lr stator self inductance, mutual inductance and rotor self inductance 

M o , M p a and p axis magneto motive force 

1 0 mj~m s moduli of three phases 

n modulus of gyrator (number of coil turns) 

n s number of effective stator coil turn 

n r number of effective rotor coil turn 

n(t) noise in time domain 

15 p p number of poles 

p probability of occurrence 

R entropy rate 

Rs, Rr stator and rotor resistances 

n-rs modulated gyrator moduh of rotor 
20 SR S , 3v stator reluctance, mutual reluctance, rotor reluctance 

5, TV average power of the signal and noise 

S« average power of the signal including noises 

s slip 

V a , Vb, V c sinusoidal input voltages 

25 x(t ) output of the degraded m?r.hm.e in time domain 

x o (/) output of the ideally healthy machine in time domain 

0) bandwidth 

<P c y <Pp a and p axis fluxes 

□ flux linkage 
30 V a „ V pj a and (3 axis stator voltages 

i^, ipj a and (3 axis stator currents 

i,^ i0 f ct and p axis rotor currents 

R,, R r stator and rotor resistances 

L,, 1^, L, stator self inductance, mutual inductance and rotor self inductance 

35 T e , Tl electro-magnetic torque and mechanical load torque 

J moment of inertia 
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c viscous resistance coefficient 

^ (0m electrical and mechanical angular velocities of the rotor 

p number of pole pairs 

X flux linkage 

5 mi-rns moduli of transformers for 3 phase to 2 phase transformation 

V a ,Vb,V b sinusoidal input voltages 

mi~m5 transformer moduli 

irt current in the rotor bar, 

m magnitude modulus that depends on the total number of bars 

I o n modulus of gyrator (number of coil turns) 

6 magnetic flux [Weber (Wb)] 

M magneto motive force [Ampere (A)] 

p permeance of circuit element [Henry (H)] 

R reluctance of circuit element [ 1 /Henry (H* 1 )] 

15 ei,e: effort 

f tl f 2 flow 

n 3 number of effective stator coil turn 

n ; number of effective rotor coil turn 

h : angular momentum [N-m-s=kg-m 2 'Sec] 

20 Rs_alpha, 
Rs_beta, 
Ri, Rsa, 

Rsb, Rsc, R electrical resistances 



25 
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APPENDIX FOR FIRST SQUIRREL CAGE INDUCTION MOTOR EXAMPLE 



In this section, for a rotor with n bars, we prove 

m= -\f ■ 

(A.1) 

5 From equations (15) and (16), the transformation matrix times its transpose is 
A T A = m 



10 



liu yivj, uoiuivimiL.uum.^ — 1 

cos* cos(* + ^) cos[e^) ... 



x m 



cos 



cos0 
cos 6 + — 

cos| 9 + — 
2(n -I)* 



0 + 



V 



sin<9 
sin 0-r — 

k "J 



(A.2) 



The result of tie multiplication is a square matrix of dimension 2, 
A T A — m 2 



S 2l S-2 



(A3) 



where 



.g^^^^j (A4) 



Equations (A.4) can be rewritten using double angle trigonometric formulas: 

i +cos <2e+^} 



15 sii- Z 
k=l 



i n-i /i_ 
"j=0 



si2=iz-{ 2 e^}=i Z 



'k=l 

n I-cos{20 + ^} 

S22= I ; 



2 j=0 

= n 1 
2 "2 



n-1 



Z cos{2e+j^} 
j=0 



(A.5) 



Via formulas 1.341-1 and 1.341-3 in Gradshteyn and Ryzhik (1980), the sums of sine and cosine terms 
on the right sides of equations (A.5) are zero, for n > 3. Thus 
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;n = l 


2 2 




l-rCOs20 


;n = 2 


K 


;n>3 


,2 




(A.6) 




i-icos2(9 


;n = l 


2 2 




■ 1-cos 2<9 


;n = 2 


n 


;n>3 


2 




(A.7) 





i sin 20 

2 

sin 28 



;n = l 

;n = 2 
;n>3 



(A.3) 

Therefore, we can conclude 



A T A = m 2 



$2\ ^ 22 



= m' 



n 


0 


2 




0 


n 






0 



■> n 
■ m~ — 



1 0 
0 1 



(A.9) 

10 for the rotor with more ttian 2 rotui bars, i.e.. a > 3. 
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Appendix For Second Example of Induction Motor 

Range of Channel Capacity Values 

In equation (27), average power of the output signal including noise is defined as 



S:=fS[x(t)ldtSj[x 0 (t) + n(tfdt 
•L c o 

* i 1 o J o 



o x o 

(A.1) 

Using equations (3) and (4), we get 
2 T 

s: = s + N'+-jfe 9 (t)-n(t)}k. 

(A.2) 

10 Shannon assumed statistical independence of *. (t) and »(0 , which made the last term of equation 
(A.2) vanish. In this article, we will remove this restriction, allowing for forms of noise 
n (/) = -K.v o (r),(0<K<l) 

(A.3) 

that can extinguish the signal, such that 

i5 *(0 = *.(0+»(0 = (i-Kk,(0 **.(0. 

(A.4) 

Here the amplitude (l - K) of x(t) diminishes with increasing noise: this can affect the signal power 
significantly. With the values just derived, 



K 



20 (A.5) 

can be less than unity, and the channel capacity can have negative, values: as noise power proportional 
K 2 increases, the output power proportional to (l -K) decreases. 



to 
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We claim: 

1 . A method of modeling a mechanical system comprising a plurality of physical components, 
comprising: 

preparing a model of a mechanical system in which at least a portion of the physical 
5 components of the mechanical system are individually modeled, wherein the model is 

configured to output data representative of the condition of the mechanical system in response 
to an input of operating conditions for the mechanical system; 

monitoring the condition of the mechanical system in response to predetermined operating 
conditions during use; 

1 0 modifying the model such that the outputtsd data of the model in response to the 

predetermined operating conditions is representative of the condition of the mechanical system 
in response to the predetermined operating conditions. 

2. The method of claim 1; further comprising predicting a failure time of the mechanical system 
using the modified model. 

15 3 . The method of claim 1 , wherein all the physical components of the mechanical system are 
individually modeled. 

4. The method of claim 1 , wherein a possible fault for each of the individually modeled physical 
components are incorporated into the model. 

5 . Thcine Lhod of claim 1 , wherein a potential failure for each of the individually modeled 
20 physical components are incorporated into the model. 

6. The method of claim 1, wherein a plurality of possible failures for the individually modeled 
physical components may interact, rendering failures not specifically associated with any single 
component, but arising from interactions between components. 

7. The method of claim 1 , wherein the condition of the modeled mechanical system is 
25 represented within die model as noise. 

S. The method of claim 1, wherein the condition of the modeled mechanical system is 
represented within the model as noise, and wherein the condition of the modeled mechanical system is 
determined by calculating a signal to noise ratio for the model. 

9. The method of claim 1, further comprising: 
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calculating the channel capacity of the modeled mechanical system, wherein the channel 
capacity is representative of the design of the system and the present condition of the 
mechanical system; 

calculating a rate of information for a predetennined job to be performed by the mechanical 
5 system, wherein the rate of information is representative of the speed, loads, complexity and 

desired accuracy of the job; and 

comparing the rate of information to the channel capacity, wherein if the rate of information is 
less than or equal to the channel capacity the model will output data indicating that the 
mechanical system is capable of performing the job at the appropriate speed, load and 
10 accuracy. 

10. A computer implemented method of modeling a mechanical system comprising a plurality of 
physical components, the method comprising: 

preparing a model of a mechanical system in which at least a portion of the physical 
components of the mechanical system are individually modeled, wherein the model is 
1 5 configured to output data representative of the condition of the mechanical system in response 

to an input of operating conditions for the mechanical system; 

monitoring the condition of the mechanical system in response to predetennined operating 
conditions during use; 

modifying the model such that the outputted data of the model in response to the 
20 predeterrnined operating conditions is representative of the condition of the mechanical system 

in response to the predetennined operating conditions. 

11. A carrier medium comprising computer instructions, wherein the program instructions are 
computer-executable to implement a method of modeling a mechanical system comprising a plurality of 
physical components, the method comprising: 

25 preparing a model of a mechanical system m which at least a portion of the physical 

components of the mechanical system are individually modeled, wherein the model is 
configured to output data representative of the condition of the mechanical system in response 
to an input of operating conditions for the mechanical system; 

monitoring the condition of the mechanical system in response to predetermined operating 
30 conditions during use; 



013270.00015:124601.01 



WO 02/27418 



PCT/US01/29935 



46 

modifying the model such that the outputted data of the model in response to the 
predetermined operating conditions is representative of the condition of the mechanical system 
in response to the predetermined operating conditions. 



5 12. A method for diagnosing a state of a system, the method comprising: 
measuring a signal from the system; 

comparing the signal from the system and an expected signal to determine a noise signal 
associated with the signal from the system; 

determining a signal strength associated with the signal from the system; 

10 detei rainin g a rate of information, the fate of rnformation associated with a desired operability 

of the system; 

detenninrng a channel capacity from the noise signal and the signal strength, the channel 
capacity being a function of a frequency spectrum of the signal; 

comparing the rate of information to the channel capacity to diagnosis the state of the system. 



15 



20 



13. The method of Claim 12 wherein the expected signal is a signal measured from an exemplary 
system operating in a known state. 



14. The method of Claim 12 wherein the expected signal is the output of a model 



15. The method of Claim 14 wherein the output of the model is adapted to approximate the measured 
signal. 



16. The method of Claim 12, the method further comprising: 
25 repeating the steps of the method over time to determine a set of diagnoses, 

detennining a prognosis of the system from the set of diagnoses. 
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